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Abstract 

The decay mode ip{2S) — ?■ J/tp + 'j'j is proposed in order to experimentally identify the effects of 
the coupling of charmonium states to the continuum DD states. To have a better understanding of 
such a two-photon decay process, in this work we restrict ourselves to investigate the contribution 
of the discrete part, in which the photons are mainly produced via the intermediate states Xcj{n-P). 
Besides calculating the resonance contributions of Xcji^P) ( J = 0, 1, 2), we also take into account 
the contributions of the higher excited states XcjC^P) and the interference effect among the IP and 
2P states. We find that the contribution of the 2P states and the interference terms to the total 
decay width is very tiny. However, for specific regions of the Dalitz plot, off the resonance peaks, 
we find that these contributions are sizable and should also be accounted for. We also provide the 
photon spectrum and study the polarization of J/V'- 
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I. INTRODUCTION 



Electromagnetic processes have always provided invaluable probes of the strong inter- 
actions, the most prominent example being the deep inelastic scattering experiments in 
the early 1970s which eventually established QCD as the fundamental theory of the strong 
interactions. Among the current challenges of QCD is the description of the plethora of 
charmoniumlike states discovered during the past years, the so-called "XYZ" (for reviews 
see, e.g., Refs. [l-4|), that lie above the DD threshold, and do not fit potential model 
expectations. It is widely accepted that the effects of the coupling of a core charmonium, 
namely, a mainly cc bound state, to a DD meson pair, the so-called coupled- channel effects, 
are an important ingredient to understand those states. We shall call the latter contin- 
uum states and the former valence states. We will argue below that two-photon transitions 
between heavy quarkonium states may provide important experimental information on the 
continuum-valence coupling. 

Historically, it has been recognized for a long time that the continuum states may shift 
the mass spectrum of a pure cc state considerably [5|-[7| . Recently, exploratory investigations 
on the mixing between discrete and continuum states in charmonium have been carried out 
in lattice QCD [8]. Lattice QCD has also provided a detailed study of the so-called "string 
breaUtag- i„ the sta* approx,mation H. wh.ch may be used to extract i,*n„at,ou on 
the continuum- valence coupling [10[. If the heavy quark mass m is much larger than the 
remaining scales in the system, then the heavy quarks move slowly in their center-of-mass 
frame, say with a typical velocity vq <^ 1, that is, with typical three momentum tuvq and 
hence with typical binding energy rnvg. Nonrelativistic QCD (NRQCD) llHl3| can be 
used to factorize the contributions from energies larger or of the order of m, and provides 
a good starting point. Recall that the heavy-light meson pair threshold lies at a typical 



nonre. 



ativistic energy ~ Aqcd, according to heavy quark effective theory (HQET) counting 



rules Then, if the binding energy of a heavy quarkonium state lies much below open 
flavor threshold, namely, ttivq ^ Aqcd one can integrate out energies of order ~ Aqcd- This 
leads to potential NRQCD (pNRQCD) in the strong coupling regime [l5|. The dynamics of 
this effective theory reduces to a heavy quark and a heavy antiquark interacting through a 



potential 



16 



- |l8j (see 191] for a review). Hence, in order to understand the properties of these 



states, there is no need to introduce exphcitly the continuum states if the potential is chosen 
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appropriately. However, for states close to or above the open flavor threshold, the coupling 
of continuum to valence states needs to be addressed, and so far it is not known how to 
proceed in a model independent way Hence, most of the analysis has been done by using 
different models. Recently, some general features of the coupled-channel effects have been 
obtained in the quark model 22| under the assumptions that valence-continuum coupling is 
described by the ^Pq model 23| and the interaction between two mesons is negligible. These 
results suggest that for the low-lying states the effect of continuum channels is hidden in the 
parameters in the potential model. In Ref . 2J| , it is mentioned that the radiative transition 



process may be sensitive to the continuum components. However, previous works, based on 



Cornell coupled-channel formalism, indicate that the relativistic corrections 25|, |26|] may be 



more important than the continuum contributions [27|]. Therefore, it is not straightforward 
to disentangle the continuum contributions in the one-photon transition process. 

Here, we propose a new process, namely, ip{2S) — )■ J/-?/' + 77 which can provide an addi- 
tional opportunity to pin down the valence-continuum coupling. From the theoretical point 
of view, electromagnetic transitions allow a cleaner analysis than hadronic decay processes. 
At the amplitude level, the contribution to this two-photon transition can be divided into 
two parts. We refer to the flrst one as the discrete part, which involves charmonium states 
only. The discrete part is dominated by the following process: the ip{2S) state decays into a 
real or virtual Xcj{nP) state by radiating one photon, and then, the real or virtual Xcj{nP) 
state decays into J/tp plus another photon. This process not only includes the cascade de- 
cay process (on-shell region), but also the off-shell region. We will study the whole phase 
space. The second part is referred to as the continuum part, in which at least one photon is 
emitted by an intermediate charmed meson. In the discrete part, the Xcj{^P){J = 0, 1, 2) 
states can be on-shell, so the total contribution of the discrete subprocess should be much 
larger than that of the continuum part. However, when the invariant mass of J/ip and one 
of the photons is far away from the resonance regions of the Xcj{^P){J = 0, 1, 2) states, the 
discrete part contribution will drop down very fast. Thus, the contribution of the continuum 
jart may be important and measurable in the off-shell region. Let us mention that in Ref. 



28|, a similar decay process of B{D)* — t- B{D) + 27 is suggested to determine the values of 



* An hadronic effective tfieory fias recently been introduced to study particular states very close to threshold, 



most notably the X(3872), which includes coupled channel effects 20 1. See also [211 for an even more 
recent proposal. 
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the strong couplings g^^^^,^^^^^ and 9^^^^,^^^^^- 

On the experimental side, such a two-photon transition process has already been studied 
sE, |30j. In recent years, more precise measurements were carried 



in the 1970s and 1980s 
out by the BESIII |3JJ 



32| and CLEO 



33 



34| Collaborations. However, they focused on the 



investigation of the cascade decay ip{2S) — )• Xcj + 7 followed by Xcj J j'^ \ ^ and on the 



study of the properties of the Xcj states. Only in Ref . |34| . was a discussion made on the 
possible amplitude of the two-photon transition in treating the backgrounds of the Xcj states. 
!Ience, so far, no one has ever used it to study coupled-channel effects. Recently, the BESIII 



35| Collaboration reported the significant data excess from the known cascade backgrounds, 
which was interpreted as the nonresonance decay of 'ip{2S) J /ip'-yy. We remark that this 
measurement is very sensitive to the line shapes of the Xcj states, especially in the data 
selection region. 

Because of the above reason, the significance of an eventual experimental determination 
of the continuous process depends very much on our knowledge on the discrete states. As far 
as we know, the discrete contribution has not been fully studied yet, and only the individual 
contribution of each XcJ state is known by using the non relativistic Breit-Wigner formula 
together with the dynamical factors to describe the Xcj hue shape in the cascade decay of 
ip{2S) — )■ XcJ + 7 and Xcj — ^ J/V' + 7 [SJ]- So, in this paper we restrict ourselves to analyze 
the discrete contribution to this decay assuming that the coupling of ip' and J /ip to DD is 
zero, and leave for a future work, the detailed evaluation of DD meson pair loops effects. 

In this work, we will use effective field theory methods to calculate the decay width, 
the photon spectrum, and the J/ip polarization in the discrete subprocess. A complete 
study of the whole contribution of the discrete Xcj{^P) states together with some higher 
radial excitations, Xcji'^P)^-, will be carried out. The rest of this paper is organized as 
foUowings: in Sec. II, we will briefly introduce the effective Lagrangian we use and determine 
the value of the effective couplings; in Sec. Ill, we will calculate the discrete contribution to 
the Tpi^S) — )► J/^/j + 77 process and show the results, and the discussion and summary will 
be given in the last section. 



The contribution of higher nP states, where n > 3, are ignored. 
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II. EFFECTIVE LAGRANGIAN FOR RADIATIVE TRANSITIONS 



The heavy quarkonium states are mainly constituted by a QQ pair and classified according 
to the spectroscopic notation , where n = 1,2,... is the radial quantum number, 

5 = 0, 1 is the total spin of the heavy quark pair, L = 0, 1, 2 . . . (or 5, P, -D . . .) is the orbital 
angular momentum, and J is the total angular momentum. They have parity P = (— 1)^^^ 
and charge conjugation C = (— 1)-^+'^. 

As mentioned in the introduction, NRQCD is a good starting point to describe this 
system. The LO NRQCD Lagrangian is invariant under S=SU{2)q(^SU{2)q spin-symmetry 
group, an approximate symmetry of the heavy quarkonium states, that is inherited in the 
subsequent effective theories. We assume that the entire dynamics of these states can be 
described by a (nonperturbative) potential. That is the case if they are in the strong coupling 



regime of pNRQCD. This is a reasonable assumption for the Xcj(l-P) states j36|, and to lesser 
extend for J/ip'^. The remaining states are close to or above threshold and are subject to 
the uncertainties due to the infiuence of DD pairs, that we plan to analyze in a separate 
work. In any case, we will only use the fact that the typical energy of the emitted photons 
is at the ultrasoft scale mt>Q, and the typical size of the system at the soft scale This 
allows hs to carry out the multipole expansion of the photon field about the center-of-mass 
coordinate, which means that the photons see the heavy quarkonium states as pointlike 
particles. Hence, it is most convenient to introduce hadronic spin-symmetry multiplets, in 
an analogous way as it was initially done in HQET jl^ . 



For heavy quarkonium states, this formalism was developed in Ref. 40|]. The states have 



the same radial number n and the same orbital momentum L can also be expressed by means 



There are indications that J /'ip may well be better described by the weak coupling regime of pNRQCD 



37 



39l |. At leading order in this regime, the dynamics is described by a (perturbative) potential, and the 



multipole expansion also holds. 
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of a single multiplet J^^^■■■^^'- [40|, 



Lv/(2L-1)(2L + 1)^^^ 
+ ■■■'^^f (1) 

where u'^ is the four velocity associated to the multiplet J'^i ■■■'^i (not to be mistaken by 
vq, the typical velocity of the heavy quark in the heavy quarkonium rest frame), K'l^'"^^ 
represents the spin-singlet effective field, and ^}J^^^■■■^^^ and ^f^^^^-^^+i represent the 

three spin-triplet effective fields with J = L — 1,L, and L + 1, respectively. The four tensors 
are all completely symmetric and traceless and satisfy the transverse condition 

v,^^-^-''^ =0 , v,^H^/-'--'' = 0, (2) 

i = 1, . . . , L, j = 1, . . . , J. The properties of H and K under parity, charge conjugation and 
heavy quark spin transformations can be easily obtained by assuming that the corresponding 
transformation rules of the multiplet J'^i - '^^- follow as: 

J^^-^"^ ^l'j,,...,,l',v^^v„ (3a) 
J''-'' A {-l)^^'C[J,,...,,fC, (3b) 

y S* J^^. ^^^S* ^, (3c) 

where C is the charge conjugation matrix (C = ^7^7° in the Dirac representation), and 
S G SU{2)q and 5" G SU{2)q correspond to the heavy quark and heavy antiquark spin- 
symmetry groups {\S,f\ = [S','^] = 0). 

Since we are going to consider the two-photon decay of ip{2S) into ip{lS) via the interme- 
diate states Xcj{nP), it will be helpful to give the explicit expressions of the S- and P-wave 
multiplets that follow from Eq.(l). For the L = S case, we have 

J='4^(H^,,-IW)K^, (4) 



and for the L = P case, 



The radiative transitions between mS and nP charmonium states in the nonrelativistic hmit 
is given by the Lagrangian: 

£ = J2^''^'"'^T^AJ{mS)J^{nP)]v^F^"' + H.c, (6) 

where 5"^'™-^ is the couphng constant, and F'^'^ is the electromagnetic tensor. The La- 
grangian in Eq.(6) preserves parity, charge conjugation, gauge invariance, and heavy quark 
and antiquark spin symmetry. 

Using the effective Lagrangian, it is straightforward to calculate the El transition decay 
widths: 

r(m S, ->nPj)- (2J + ^3^^^^ , (7a) 

linPj^mb^)- ^^^^3^ , (7b) 

le energy of the emitted photon. The results are slightly different from those 



where k-y is t 

given in Refs.|40l. l4l|. It is because the initial and final charmonium states can not be static 
simultaneously, and we choose different values of for them to maintain the transverse 
condition in Eq.(|2]). Namely, in the vertex ([6]) we have substituted the velocity v'^ in the 
current that produces the outgoing particle by Vp its four velocity. The remaining velocities 
(i.e. the explicit one and the one in the current that annihilates the incoming particle) are 
chosen as ff , the four velocity of the incoming particle. The explicit calculations in ([7]) 
show that these changes maintain the spin-symmetry ratios between the decay widths of 
the different states ( provided that the same mass is used for each spin multiplet ). Note 



that Vi.Vf = 1 + 0{vq), and the results of Ref. j4l| differ from ours by 0{vq) contributions. 



In ([7]) we have also allowed the coupling constant S"^'^^ of the spin-symmetry multiplet to 
depend on the total angular momentum J, namely, S"^^'"^^ — S^^'"^^ = §nP,mS _j_ (9(1/777,2^ _ 
In this way we are accounting for spin-symmetry breaking terms due to the spin-orbit, 
spin-spin and tensor potentials in pNRQCD (or in potential models). Note that although 
spin-breaking terms look like an 0{l/m) in NRQCD (like in HQET), this is actually the 
case only when ultrasoft particles are emitted, for instance, in magnetic transitions, or in 
processes involving pseudo-Goldstone bosons [21]. Soft emissions can only be virtual and 



lead to the above mentioned 0(l/m^) spin-breaking potentials, which make the coupling 
constants different at this order. In conventional potential models, the differences between 



rnP.mS 



are usually taken to be of order fg, although on general grounds they could be 



as large as order vq [17|. ^'^P''^^ can, in principle, be obtained by calculating the matrix 
element of the electromagnetic current between the wave functions of the nP and mS states 
in pNRQCD (or in any potential model, see Ref. 42(| for a recent review). Nowadays they 



can also be obtained from lattice QCD 



43 



44| . Here we will determine their values using 



experimental data, except for those related to the Xcj(2P) for which experimental data is 
not available. In the last case the values will be estimated using a potential model. 

For the IP case, the '>p{2S) — > Xcji)-P) + 7 and Xcj(l-P) J j'^ processes have been 
measured with a very high precision [45], so we can determine the values of the relevant 



coupling constants by comparing with the experimental results. Using the PDG data 



25 



V(2S) 



3.686GeV, M 



V.(15) 



3.097GeV, M. 



XcJ 



3.415,3.511,3.556 GeV (for J = 0,1,2), 



304keV, r^^^ = 10.3, 0.86, 1.97 MeV (for J = 0, 1, 2), and 

Br(^(2^) ^ Xco + 7) = 9.62% ; Br(xco ^ Jl^ + 7) = 1-16%; 

Br(^(25) ^ Xci + 7) = 9.2% ; Br(xci ^ J/^ + 7) = 34.4%; 
Br(V^(2S) ^ Xc2 + 7) = 8.74% ; Br(xc2 ^ W + 7) = 19-5%; 



(8a) 

(8b) 
(8c) 



we obtain the absolute values of 5^^'^"^ and b^F'^^ listed in Table I. For the IP case, only the 



Xc2(2P) [formerly called Z(3930) 46|,] has been included in the PDG, so we estimate the 
relevant parameters with the help of a potential model. The spectrum of the 2P states and 
their radiative decay into the lower S-wave states have been calculated by numerous groups 
(for reviews see, e.g., Ref.[42|]). By using the screened potential, the updated results of the 



charmonium spectrum and the El transition rates are given in Ref. 26|] . In this work, besides 
setting M^^2(2p) = 3.929GeV, X(3872) is assigned to the Xci(2-P) state and M^^„(2p) = 
3.842GeV is predicted. From the nonrelativistic results of r(xcj(2P) — ?■ J /'ip{ip{2S)) + 7) 
presented in Table IV of Ref. ^6| , we obtain the absolute values of the corresponding coupling 
constants and give them in Table I. Because the mass differences between the 2P states are 
much smaller than the mass gap between the 2P and 25* (15*) states, in practice we use the 
center of gravity mass 

M^,o(2P) + 3 X M^^^(2P) + 5 X M^^,(2P) 



M, 



(2P) 



9 



(9) 



8 



TABLE I: The numerical values of the coupling constants 6j^'"^^{GeV~^) are shown. For the 
n = 1 case, the results are obtained by fitting the experimental data, and for n = 2, the results are 
determined by comparing with the potential model predictions 26|. 





Xco{lP) 


Xcl{lP) 


Xc2{lP) 


Xco{2P) 


Xci{2P) 


Xc2(2P) 




0.211 


0.230 


0.228 


5.27 X 10-2 


5.30 X 10-2 


5.34 X 10-2 




0.224 


0.235 


0.273 


0.410 


0.413 


0.416 



and the center of gravity coupling constant 

!:2P,ms _ "O 



8' 



r2P,mS I Q ^2P,mS , r- ^2P,mS 
On + 6 Oi + 5 X ^2 



(10) 



rather than calculating the contribution of the individual 2P states. 



III. DISCRETE CONTRIBUTION TO ^(25) ^ J/ip + 27 

Now, we proceed to calculate the decay rate of the process '?/'(2S')(po) — >■ •^/'^(Pi) +7(^*2) + 
7(^3) via intermediate states Xcj{nP). Such 1 — 3 process can be described by the following 
dimensionless variables: 

2po ■ Pi 



Xi 



E^^ = 2' (11) 



In terms of Xi, the three-body phase space $(3) can be written as 

^l(2S) 

d$(3) = ^^^^^3 (5(2 -X1-X2- X3)dxidx2dx3, (12) 

For each intermediate Xcji'^P), there are two Feynman diagrams, which are shown in Fig.l. 
The corresponding Feynman amplitude is denoted by Putting the contributions 

of the three 2P states together, the total Feynman amplitude is then divided into four parts: 

Each of Ai on the right-hand side of the above equation can be obtained from the La- 
grangian in Eq.(6) upon making the same replacements as discussed after ([6]). Since the 
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liPi) liPa) 7(P3) 7(P2) 






Xr.iinP) Xrj(nP) 
FIG. 1: The Feynman diagrams for 'ip{2S) decay into J/ip + 27 via intermediate states XcjinP). 

propagator of the IP fields may become on-sliell, self-energy corrections must be consid- 
ered. We approximate them by introducing a constant decay width, which, parametrically, 
is 0{ma'^VQ, mavg). Being that these figures are much smaller than mvg, 1/m corrections 
to the static propagator must be considered in order to match its size. As an alternative, 
we use relativistic propagators, that include them. For example, the Feynman amplitude of 
the first diagram in FIG. [T] for Ai^'''^^^^^ which is the most complicated one, is 

^ ^ 2 25 2 V'^-I + I* r;^,2(iP)/M^,2(lP) 

X Tr[i±4^«.^;^i±^] X v^,2sv.,,FrFr. (14) 

where W^^ = {-g^^ + v^), Pf = p^ef - pf e^, = {p'o - P^)/M^..(ip), and r^^,(ip) is 
the total width of Xc2{^P)- In the above expression, we have omitted the imaginary unit 
I, which is a global factor and has no influence on the final result. For the convenience of 
further discussion, we also divided the decay width of the discrete part into four parts: 

T,u^{2S) ^ j/^ + 77) = tZ + rZ + ± r^^^^^ (is) 

where is the sum of the three individual contributions of the Xcj states, which 

is proportional to X]j=o l-^^''"'''^^^!^' -'^int the interference between the IP states (~ 
M*^''-^''^^^M^'=J'^^^^), F^^ is the contribution involving the 2P states only, and 
Fj^^^ is the interference contribution between the IP and 2P states. The "+(—)" sign before 
the last term corresponds to the two possible relative phase angles (vr) between the IP 
and 2P states. 

We have computed ^ |A^™p analytically, but the outcome is too lengthy to be presented 
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— Xj'^P) mode 
---- Xj'^P) mode 




CM 

X 









lnt(xJ1P),xJ1P)) 

Int(x (1P),xJ1P)) 

mxJ^P),xJ^P)) 




■ 1 

• /V. ■ 


: 1 1 


1 .... 1 ... , 



0.1 



0.2 



0.3 



(a) 



(b) 



FIG. 2: The partial decay width as a function of the photon energy fraction X2- (a) the individual 
contribution of the three Xcj {J=0,l,2) states, corresponding to Fj^^ in (fT5]) . (b) the contribution 
of the interference terms between the three Xcj (<^=0,1,2) states, corresponding to Fj^^ in (fT5]l . 



here. After doing the phase space integrals numerically, we obtain that 

rZ = 15.14keV ^ Yl r(^(2^) ^ 7 + Xcj) x Bt{xcJ ^ W + 7), 

4.13 X 10" VV. 



r^P K OK XX in-2u„-\A ■p2P o on XX in— 3i,„-\a -pliSf /i i q xx i n-2i 

Int 



5.95 X lO^^keV, F^^ = 2.80 x lO^^eV, Fj^^ 



(16) 



The numerical results yield that on the total decay width level the effects of the interference 
among the Xcj(l-P) states as well as the effect of the 2P states are so small that they can be 
neglected. We have also calculated the photon spectrum d$(3)/(ix2, and display the figures 
for each part separately in FIG. [2] and FIG.[3l respectively. 

Since we have chosen different v for the initial and final states in Eq. (fT^ . we may also 
have a nontrivial J/ ip polarization, which should be zero in the single v case. Similar to the 
production case 47|, here, we define the polarization parameter a as 



a 



(17) 



rT + 2Fi' 

where F^ and F^ are the decay widths for J/ ip in the transverse and longitudinal polarization 
states, respectively. We calculate a in the rest frame of ip{2S) and get a ^ —0.16. This is 
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0.1 



0.2 



0.3 



(a) (b) 

FIG. 3: The partial decay width as a function of the photon energy fraction X2: (a) the contribution 
of the IP states, corresponding to in (jlSp . (b) the contribution of the interference terms between 



the 2P and the three IP states, corresponding to ^{^^ in (fl5D . 



slightly different from zero, and we interpret it as a purely kinematic effect. We find that a 
is mainly determined by the resonance contribution, being the influence of the interference 
and the 2P terms tiny. 

The experimental measurement of the 'ip{2S) — ?■ J/ip + 'yy decay can be implemented in 
the off-mass shell region of Xcj, i-e., the experimentally sensitive region in Daltiz plot 

0.15 < M^^ < 0.51 GeV, 3.43 < Mj/^^ < 3.49 GeV, (18) 

where M^^ is the invariant mass of the two photons and Mj/^^ is the invariant mass of 
J/ip and the higher energy photon. These cut can mostly exclude the contribution from 
the highly yielded Xcj{^P) states. From the photon spectrum in the cut region indicated in 
FIG. m the different contributions to the decay width ( |T5|) read 

Fi^fd = 4.68 X lO^^eV, Fj^^ = 6.5 x lO^^keV 
F^^ = 1.82 X 10~^keV, Fj;^^ = 4.78 x lO^^keV. (19) 

We have also computed the branching ratio, photon spectrum, and the polarization param- 
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I I I I I I I I I I I I I I I I nlTT-r-r~~i— i~i~i~TTTT~r~;~'i ni"' ■ 

0.11 0.12 0.13 0.11 0.12 0.13 0.11 0.12 0.13 

X2 Xj X2 

(a) (b) (c) 

FIG. 4: The discrete contributions to the photon energy spectrum of the ip{2S) — )• J/tp + 77 
process in the cut region: (a) the contribution of the IP states, corresponding to Fj^^ and Fj^^ in 
(jl5p . (b) the contribution of the 2P states and of the interference terms between IP and 2P states, 
corresponding to F^^ and to Fj^^'^ in ()15p . (c) the total contribution for a different relative phase 
angle 0, corresponding to the it sign in (|15p . 



eter in the cut region. The result of the branching ratio is 



BX{ij{2S)^J/^ + 2^) 



1.92 X 10~^ (for ^ = 0) 
1.60 X 10-^ (for e = n) 



(20) 



The result in Eqs. fll9ll20p shows that, contrary to the case of the total decay width (fT6l) . in 
the cut region the effect of interference among IP states as well as the contribution from the 
2P states can not be ignored. It is more than 10% of the sum of the separated contributions 
of the Xcj(l-P) states. The J/ip produced in the cut region tend to be unpolarized and the 
polarization parameter in the cut region becomes a^^* = —0.122 and —0.107 for 6' = and 
6 = IT, respectively. If we only include the three individual contributions of Xcj(l-P), the 
value of a'^"* turns to be —0.078. Finally, the different contributions to the photon spectrum 
in the cut region are shown FIG. HI 

As mentioned before, the coupling constant S"^^'"^^ is related to the spatial matrix element 
{nP\r\mS) = Rnp{r)Rmsij')^^dr in pot ential models. At least three different potential 



models, including the Cornell potential and the screened potential [26|, give the phase 



angle ^ to be vr 



AW\. Hence, the 6 = n option in our calculation appears to be favored. 
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IV. DISCUSSION AND SUMMARY 



We have estimated the discrete contribution to the il'(2S) J/V' + 77 due to electric 
dipole transitions in the whole phase space and, in particular, in the cut region used at 
BESIII. Higher multipole electric transitions and the magnetic ones are suppressed by at 
least Vq in the amplitude. This is also the case for contributions arising from two-photon 
vertices. The largest uncertainty in our calculation comes from the fact that we neglect the 
contribution of nP states for n > 3. This contribution is not parametrically suppressed 
by powers of vq, although we expect it to be small for the following two reasons: (i) the 
propagator of the nP state is increasingly off-shell when n increases, and (ii) the coupling 
constants 6^^'^^ and Sj^'"^^ are proportional to the overlap of the radial wave functions of 
the corresponding states, which decreases with n. Apart from this truncation, the estimate 
is reliable at leading order in Vq. In fact, the modifications we have made in (fT^ . which 
introduce terms of higher order in vq, take into account relativistic effects in the kinematics, 
and hence help in providing a better estimate. 

From the point of view of potential models, the leading corrections to the El transition 
processes, 'ipi^S) — )■ Xcj+1 and Xcj J/i^+l, mainly arise from the three following sources: 
(i) relativistic modification of the nonrelativistic wave functions, (ii) finite size effects, and 
(iii) contribution of high Vq order electromagnetic operators [42,] . In this work, we use the 
effective Lagrangian, Eq. (fT^ . to describe the El transition process. For the IP case, we 
determine the values of the corresponding coupling constants from the experimental data. 
Thus, the effects of (i), which are dominant, are taken into account in our results, as far as 
the vertices involving the IP states is concerned. Therefore, our results for the two-photon 
decay width in the off-shell region are more accurate than a potential model calculation for 
the one-photon El transition process. Both the corrections due to (ii) and to (iii) could 
be taken into account by including higher dimensional operators in the vertices (|6]). As 
mentioned before these operators are suppressed by at least order Vq. 

Let us next discuss how our results compare with the usual inputs in the Monte-Carlo 
(MC) codes that are used to analyze the experimental data. In the experimental treatment 
of the four contributions in Eq.( [T5|) in the peaking region, usually only the first one, Fj^^, is 
taken into account, and is often modeled using the nonrelativistic Breit-Wigner line shape 
of Xcj and J/ip. The other three components are negligible and generally omitted. But in 
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the off-shell region, as argued above, the other three components will be sizable and have 
to be considered in the data treatment. Furthermore, in the off-shell cut region f|T8|) . even 
this naive nonrelativistic Breit-Wigner line shape description of the individual Xcj states 
contribution needs improvement. The reason is that a (single) nonrelativistic Breit-Wigner 
is only a good approximation to the line shape of dr^^/dx2 in the resonance peaking region. 



as shown in FIG. 2(a) In the cut region equation (fTSl) . which lies between the Xco and Xci 
resonance peaks, there is no guarantee that the nonrelativistic Breit-Wigners will provide a 
correct description of data. Let us next point out the main ingredients it misses. In the El 
transition process, the decay rate is proportional to the factor as shown in Eq.([7]). We 
may then improve on the nonrelativistic Breit-Wigner approximation by just introducing 
the correct photon energy dependence in each vertex i.e. the full k^ik^2 scale factor from 
the two El transitions. Note that if one only includes the k^ correction due to the first 
El transition, F^"^ in the cut region will be overestimated, because the energies of the two 
photons are negatively correlated. 

In FIG. [5l the MC simulation results of the line shapes in the cut region are shown, which 
are implemented in two ways. One is done by including the k^^k^2 correction from a double 
El transition, and the other one is not. The difference between these two MC simulation 
results can be understood as follows: the El transition enhances the right tail of photon 
energy peak and depresses the left tail of the nonrelativistic Breit-Wigner. As a whole, 
the effect of the correlated emitted photons in the double El transition increases the XcJ 
contribution in the cut region. For comparison, we also plot the effective Lagrangian result 
of the sum of the three individual XcJ states contribution in FIG. El It is clear that neither 
of the simulation results agree with that of the effective Lagrangian calculation, although 
they are qualitatively similar. From the amplitude in Eq. ffT^ . it can be verified that the 
^71^72 factor and the nonrelativistic Breit-Wigner is only the leading-order nonrelativistic 
approximation of the effective Lagrangian calculation in the off-shell region. Hence, even 
for a delicate description of the individual contribution in the off-shell region, including the 
double El transition correction in the naive nonrelativistic Breit-Wigner only may not be 
enough. 

In summary, we have estimated the discrete contribution to the ip{2S) J/ip + 'J'J due 
to electric dipole transitions in the whole phase space, and in particular in the cut region 
in the experimental measurement. We find that for the full decay width the interference 
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FIG. 5: The MC simulation of the cascade decay of ip{2S) — )• (</''/'7i)xcj72 ™ the cut region, where 



the branching fractions are from PDG 



45j . The dotted hne denotes the naive nonrelativistic Breit- 



Wigner simulation, the solid line is the simulation including the k^^ k^^ factor, and the dashed line 
is the contribution of the three individual Xcj(l-P) states, corresponding to Fj^^^ in (fTSl) . calculated 
in this paper. 



contribution and the contributions of higher excited states can be safely neglected. However, 
in the regions of the phase space off the resonance peaks, their contributions are considerable 
and important for a delicate experimental measurement. As argued in the Introduction, a 
large deviation of our results from an experimental observation, would indicate that the 
effects of the DD threshold are significant. 
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